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■ We calculate the distribution of the current density j in superconducting films along the direction 

of an external field applied perpendicular to the film plane. Our analysis reveals that in the presence 
of bulk pinning j is inhomogeneous on a length scale of order the inter vortex distance. This 
inhomogeneity is significantly enhanced in the presence of surface pinning. We introduce new critical 
state model, which takes into account the current density variations throughout the film thickness, 
P3 . and show how these variations give rise to the experimentally observed thickness dependence of j 

O ' and magnetic relaxation rate. 
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I. INTRODUCTION 



The magnetic behavior of type-II superconductors depends strongly on the sample shape . Significant progress 
'■^ , has recently been made in understanding the effects of a sample aspect ratio on its magnetic behavior jl], in particular, 
^ ' in the case of thin films with the magnetic field normal to the film plane ( "perpendicular geometry" ) |^,|| . Theory 
O a-iid experiment show that the magnetic behavior in the perpendicular geometry has many distinctive 

features, essentially different from the parallel geometry, e.g. a more complicated structure of the critical state and 
' the presence of geometrical barriers . 

A number of elegant analytical solutions for the perpendicular geometry (for strips and disks) describe the Meissner 
^ 1^, the mixed state and magnetic flux creep These solutions arc based on the important ansatz, that 

one can treat the film as an infinitesimal thin plane. Then, current distribution related to vortex bending does not 
infiuence the results of the analysis which deals only with the current density and the vortex displacements averaged 
, over the film thickness. This approach was very successful in explaining the peculiarities of the current density and 
• the magnetic induction distribution across the film plane. However, this approach cannot account for any thickness 
00 dependence of both persistent current density j |p|,p|-p^ and magnetic relaxation rate ^1 12 1 in thin films. 



On 1 Explanation of the observed decrease of j with the increase of the film thickness d is usually based on the idea that 
pinning on surfaces perpendicular to the direction of vortices is strong enough and must be taken into account • 
However, as we demonstrate below, this is not sufficient for understanding the thickness dependence of the magnetic 
relaxation rate, which was found to decrease with the increase of the film thickness pd] , p^ . 

Another explanation of the observed thickness dependence of the current density may be based on collective pinning 
(— ( in a 2D regime, i. e., for longitudinal correlation length L larger than the film thickness. This case is carefully 
Q considered in | |T^ . In a 2D collective pinning regime, the pinning is stronger for thinner samples. As a result, in this 
O ■ model both the critical-current density and the creep barrier are larger in thinner samples, contrary to the experimental 
] results. Also, this scenario is probably not relevant for the explanation of the experimental data discussed below, 
because the thickness of our films d > 800A is larger than L « 40 - lOOA. 

In order to understand the experimental results we calculate the current density and magnetic induction distribution 
by using the 'two-mode electrodynamics' theory suggested earlier to explain the AC response in bulk materials 
. The essence of this theory is that two length scales govern the penetration of fields and currents into type-II 
superconductors. The longer scale is of electrodynamic origin and, therefore, is more universal: it exists, for example, 
in a superconductor in the Meissner state (the London penetration depth) or, in a normal conductor (the skin depth). 
The shorter scale is related to the vortex-line tension, so it is unique for a type-II superconductor in the mixed state. 
This scale was introduced into the continuous theory of type-II superconductors by Matheiu and Simon [|l6| (see also 
| p7| , p^ ). When applying the two- mode electrodynamics to the critical state one may ignore the time variation, i.e., 
the two-mode electrodynamics becomes the two-mode electrostatics theory. 
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Our analysis of a type-II thin superconducting film within the two-mode electrostatics theory leads to the conclusion 
that for strong enough bulk pinning, inhomogeneity of the current density becomes important, even in the absence of 
surface pinning, if the film thickness exceeds the Campbell penetration depth Ac- Thus, inhomogeneity of the current 
distribution throughout the film thickness is a distinctive and inevitable feature of the perpendicular film geometry 
like, for example, the geometrical barrier j^]. Inhomogeneity of the current distribution is significantly enhanced if 
the critical state is supported by the surface pinning. In this case, most of the current is confined to a layer of a 
depth of the order of the intervortex distance, which is usually much smaller than the London penetration depth A 
and film thickness. As a result of this inhomogeneity, the measured average critical current density becomes thickness 
dependent. This current inhomogeneity also causes a thickness dependence of the magnetic relaxation rate. In the 
following we present detailed calculations of the distribution of the current density j and induction field B in thin 
type-II superconducting film, resulting from surface and/or bulk pinning. We then introduce the first critical state 
model which takes into account the variation in j throughout the film thickness. Calculations based on this critical 
state model lead to a thickness dependence in j and magnetic relaxation rate. These predictions are compared with 
the experimental data. 



II. THEORY 



A. Equations of electrodynamics for the mixed state in perpendicular geometry 

Let us consider a thin superconducting strip, infinitely long in the y -direction, with width 2w {—w < x < w) and 
thickness 2d (— d < z < d). External magnetic field H is applied along the z -axis, perpendicular to the film plane. 
The vortex density n is determined by the z-component of the average magnetic field (magnetic induction) B 
in the film: n = B^/^o- Supercurrent of density ly {x,z) flows along the y-axis resulting in a Lorenz force in the 
x-direction, and a vortex displacement u along the x-axis. 

We begin with the electrodynamic equations describing the mixed state of type-II superconductors in such a 
geometry. They include the London equation for the x-component of the magnetic field: 

^d^B^ _ du 



the Maxwell equation: 



and the equation of vortex motion: 



where 



47r , _ dB^ dB, 



"^dt c^''^ An dz^ ' 



H* = ^ln^ (4) 



is a field of order of the first critical field ifci, ao — y^o/Bz is the inter- vortex distance, and rc ~ ^ is an effective 
vortex core radius. The equation of the vortex motion arises from the balance among four terms: (i) the friction force 
proportional to the friction coefficient 77; (ii) the homogeneous, linear elastic pinning force ex k (i. e. assuming small 
displacements u); (iii) the Lorentz force proportional to the current density j; and (iv) the vortex- line tension force 
(the last term on the right-hand side of Eq. (^)), taken from Ref. |T^ ]. 

In the parallel geometry, {d — > 00), vortices move without bending so that the x-component Bx is absent, and 
the Maxwell equation becomes: A-Kjy/c — —dBz/dx. Since B^ is proportional to the vortex density, this current 
may be called a diffusion current. The case of the perpendicular geometry, [d <C w), is essentially different: the 
diffusion current is small compared to the bending current dBx/dz (see the estimation below) and may be neglected 
for calculation of the distribution throughout the film thickness (along the z-axis). As a result, Eq. (^) becomes 

'^at 4^ dz +4^^ dz^ ■ 



2 



Equations (p and (||) determine the distribution of the displacement u{z) and of the in-plane magnetic induction 
Bx{z). This also yields a distribution of the current density {AiT/c)jy{z) = dBx{z)/dz. But these equations are still 
not closed, since the two components of the magnetic induction, and Bz, and current density jy{z) are connected 
by the Biot-Savart law. However, neglecting the diffusion current in the Maxwell equation we separate the problem 
into two parts: (1) determination of the distribution of fields and currents along the z- axis, taking the total current 
ly =^ cB^/2-K (here B^. = B^ (z = d)) and the perpendicular magnetic- induction component B^ as free parameters; 
(2) determination of the parameters ly and Bz using the Biot- Savart law. The latter part of the problem (solution of 
the integral equation given by the Biot-Savart law) has already been studied carefully in previous works . In the 
present work we concentrate on the analysis of the distribution of fields and currents throughout the film thickness 
(^-dependence). 

The accuracy of our approach is determined by the ratio of the diffusion current dB^/dx to the bending current 
dBjc/dz, since we neglect the diffusion current contribution to the total current. Suppose, as a rough estimation, 
that Bz ~ Bx 1^. Then, the diffusion current density is roughly ~ ly/w, whereas the bending current density is 
~ ly/d Thus, the ratio between the diffusion and the bending current is approximately d/w ~ 10^'^ lO^'^ 

for typical thin films. Note that this condition does not depend on the magnitude of the critical current and is well 
satisfied also in typical single crystals, where d/w ^ 0.01 0.1. Therefore, the results we obtain below hold for a wide 
range of typical samples used in the experiment. 



B. Two-mode electrostatics: Two length scales 

Let us consider the static case when vortices do not move, hence there is no friction. Then, Eq. (^ becomes 

Ait oz Ait oz'' 

Excluding the B^ component of the magnetic induction from Eqs. (|^) and (H) we obtain equation for the vortex 
displacement: 

The two length scales which govern distributions over the 2:-axis become evident if one tries to find a general solution 
of equation ^ in the form ^ u ^ exp(ipz). Then, the dispersion equation for p is bi-quadratic and yields two 
negative values for p^. In the limit k ^ AttX^ /^q{H* + Bz) (weak bulk pinning): 

2 1 1 H*+Bz 

2 1 47rfc 

^'^'^^'MH* + Bz) ' ^ ' 

Thus, the distribution along the 2— axis is characterized by the two length scales; the Campbell length Ac, which 
is the electrodynamic length, and length A, given by Eq. (H), which is related to A and the vortex- line tension. 



C. Current density and field distribution 

In order to determine distribution of currents and fields throughout the film thickness, one must add the proper 
boundary conditions to the general solution of Eq. (|^). We look for a solution which is a superposition of two modes. 
In particular, for the vortex displacement we can write: 

z z 

u{z) — Uq cosh 1- ui cosh — . (10) 

A 

Using Eq. (^) one has for the current density: 

477. dB^ Uq z Ui z 

-'y-^^ A^ "^"^^ " A^ A • ^ ^ 
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The total current is 

—ly = 2BJd) = 2B,^ sinh - 2i/*^ sinh £ . (12) 
c Ac Ac A A 

Equation (|l2|) is in fact a boundary condition imposed on the amphtudcs of two modes, uq and ui. The second 
boundary condition is determined by the strength of the surface pinning. If displacements are small, the general form 
of this boundary condition is 

= , (13) 

±d 

where a = in the absence of surface pinning and a — )■ cxd in the limit of strong surface pinning. In the following 
parts of the section we consider these two limits. 



au(±d) ± 



du 



1. Surface pinning 

Let us consider the case of surface pinning in the absence of bulk pinning [k — 0), when the Campbell length 
Ac — > oo (see Eq. (^). By "surface pinning" we understand pinning due to surface roughness on the surfaces 
perpendicular to the vortex direction. The surface roughness is assumed to be much smaller than the film thickness 
d. By substituting Ac — s- oo in the general solution Eq. (Eoh, we derive the displacement for surface pinning: 



2; 

u{z) = Uq +ui cosh — , (14) 
A 



where uq and ui are constants, which can be determined from the boundary conditions Eqs. (12) and (|l3|). Note, 
however, that uq is not important in the case of surface pinning, because the constant uq does not affect distributions 
of currents and fields. 

The magnetic field Bx is obtained from Eq. (^): 

Bjz) = -H*^smh^, (15) 
A A 

and the current is determined from the Maxwell equation (^) neglecting the diffusion current: 

j„ = -^ff*^cosh^. (16) 

It is important to note that the characteristic length A , which varies between the London penetration length A 
and the inter- vortex distance ao ^ y^^o/Bz, is much smaller than A for a dense vortex array, B^ 3> H* . Taking 
into account that usually thin films have thickness less or equal to 2A, the effect of the vortex bending due to surface 
pinning may be very important: most of the current is confined to a thin surface layer of width A. 

The current density on the surface is js = jy (z = d) = — (c/Air) H* ^ui/A^^ cosh ^d/ A^ . Thus, 

c H* cosh ■ 



"1=-- TTTT-^- (17) 



A 

The total current integrated over the film thickness 2d is: 

ly = [ jy{z)dz — — sinh — = 2Ajs tanh — . (18) 

J-d 27r A A A 

Thus, the average current density ja = Iy/2d - the quantity derived in the experiment - decreases with thickness as 

A , d , , 

Ja=Jsjtanh^, (19) 
d A 



yielding ja = jsA/rf for X/d << 1 as found experimentally 
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The field and the current distribution over the film thickness are: 



T cosh f cosh f 

J Jz) = ^ = J, ^, (20) 

•^"^ ' 2Asinh| ^ cosh I' ^ ' 

9_- sinh 5. zItt ~ sinh i. 
c smh S c coshS 

A A 

Thus, the current penetrates into a small depth A and is exponentially small in the bulk beyond this length. 

2. Bulk pinning 

A remarkable feature of the perpendicular geometry is that, even in the absence of surface pinning, vortices are 
bent. This is in striking contrast with the parallel geometry where the diffusion current distribution is homogeneous 
along the direction of vortices and, therefore, does not bend them. Absence of surface pinning means that at the 
surface du/dz = {a, vortex is perpendicular to an ideal surface). This yields the relation between uq and ui [see Eq. 

A sinh ^ 

Ul = -uo- . ■ ^ 

Ac smh J 



Then, Eq. (pa) becomes 



i^/, = 2(S. + ff*)^sinh-^ 
c Ac Ac 



ly = 2{B^ + H*)^ sinh — . (22) 



The current distribution is 



/ 1 cosh ^ 1 H* cosh | \ 

Jyi^) = ly [jx^,H* + B, sinh 1 + 2XH*+B^smh i ) " ^^^^ 



In the limit d <C Ac Eq. (E3) yields 



IB, 1 H* cosh.^ 



Jyi^) -^yy2dH*+B^+ 2A H* + B^ sinh f ) ' ^^^^ 

Another interesting case is that of the dense vortex array, B^ ^ H* : 

L cosh cosh 

Jv iz) = ^ ^ = js ¥ , (25) 

' 2Acsinh ^ cosh^ ^ ^ 

Ac Ac 

where again is the current density on the film surface. Remarkably, current density is inhomogeneous even in 
the absence of surface pinning. We illustrate this in Fig.l, where we plot jy (z) /jf, vs. z/d at different ratios d/Xc- 
"Uniform" bulk current density jb — Iy/2d corresponds to the limit d/Xc — 0. Physically, such current profiles reflect 
Meissner screening of the in-plane component B^ of the self-field. 
For the average current density we have 

Ac d 

ja = js —r tanh — , (26) 
a Ac 



which is similar to the case of the surface pinning, Eq. (^^, with A replaced by Ac. 

Thus, in the perpendicular geometry, the current distribution is strongly inhomogeneous: the whole current is 
confined to a narrow surface layer of width A (surface pinning) , or Ac (bulk pinning) . 
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D. Critical state 



In the theory given in the previous sections we have assumed that currents and vortex displacements are small. In 
this section we deal with the critical state when the current density equals its critical value jc- Let us consider how 
it can affect our picture, derived in the previous sections for small currents. 

1. Surface pinning 

If vortices are pinned only at the surface, the value of the critical current depends on the profile of the surface, 
and one may not use the linear boundary condition imposed on the vortex displacement, Eq. ( |l3| ). However, the 
^-independent vorte x displacement uq does not influence the current density and field distribution in the bulk as 



shown in Sec. II C 1 (see Eqs. (g_5|) and (|16|)). Therefore the bulk current density and field distribution derived from 



our linear analysis can be used even for the critical state. 

2. Bulk pinning 

In this case our theory must be modified for the critical state. In particular, for large currents the bulk pinning 
force becomes nonlinear and, as a result, the current and field penetration is not described by simple exponential 
modes. Formally, this nonlinearity may be incorporated into our theory assuming a u - dependent pinning constant 
k, thus allowing k to vary along the vortex line. As an example, let us consider the case of strongly localized pinning 
force when the vortex is pinned by a potential well of a small radius like that sketched in Fig. 2: the vortex energy 
per unit length (vortex-line tension) is given by e for vortex line segments outside the potential well and by £o for 
segments inside the well. Thus, the pinning energy per unit length is e — Eq- In fact, such a potential well model may 
describe pinning of vortices by, for example, one-dimensional columnar defects or planar defects, such as twin or grain 
boundaries |l^j2^. The latter is relevant in thin films obtained by usual method of laser ablation. Therefore, we can 
also use such a pinning potential as a rough qualitative model for typical types of pinning sites, in order to illustrate 
the effect of bulk pinning on the current density distribution and the rate of magnetic relaxation in thin films. 

If the current distribution were uniform, such a potential well would keep the vortex pinned until the current density 
jy exceeds the critical value c(e — eo)/$o''d- The escape of the trapped vortex line from the potential well occurs via 
formation of the un-trapped circular segment of the vortex line (see Fig. 3(a)). In this case, both the critical-current 
density and the energy barrier for vortex depinning do not depend on film thickness |]l9| . 

But, in perpendicular geometry the current distribution is not homogeneous. In order to find it for the critical 
state, we may use the following approach. The vortex line consists of the trapped and untrapped segments as shown 
in Fig. 3(b). The untrapped segment is beyond the potential well, therefore there is no bulk pinning force act ing on 



it. This means that the shape of this segment is described by Eq. (||) with fc = 0. Applying the theory of Sec. II C 1 



one obtains that the total current ly = J^j^jy{z)dz is concentrated near the film surfaces within a narrow surface 

layer of width A. Inside the surface layer the vortex line is curved, but has a straight segment of length L outside the 
layer, as illustrated in Fig. 3(b). As for the vortex-line segment trapped by the potential well, we assume that it is 
straight and vertical, neglecting its possible displacements inside the potential well. Formally speaking, our approach 
introduces a non-homogeneous bulk-pinning constant k assuming that fc = for the untrapped segment and k — oo 
for the trapped one. The energy of the vortex line in this state is determined by the line tensions (e and Eq) and is 
given by 

L $0 f ^0 \ 

E = 2e 2eoL ~ 2 — /uLtano; = 2Ltana esina ly ] , (27) 

cos a c ^ \ c J 

where the contact angle a is determined by the balance of the line-tension forces at the point where the vortex line 
meets the line defect: 

cos a = — . (28) 

£ 
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III. MAGNETIC RELAXATION 



We now discuss the effect of current density distribution on the thickness dependence of magnetic relaxation. We 
first show below, that uniform current density cannot explain the experimentally observed thickness dependence. We 
also show that inhomogeneous current density distribution, resulting from the surface pinning only, cannot explain the 
experimental data too. We demonstrate that only presence of a bulk pinning and the resulting current inhomogeneity 
may lead to an accelerated relaxation in thinner films. We finally discuss the general case when both bulk and surface 
pinning are present. 

As pointed out above, if the current distribution is uniform throughout the film thickness, a trapped vortex may 
escape from the potential well (Fig. 2) via formation of a circular segment of the vortex line (Fig. 3(a)), with the energy 

E = eL- eoLo - ^jyS, (29) 

where L and Lq are the lengths of the vortex line segment before and after formation of the loop, S is the area of the 
loop pO| , p^ . If the loop is a circular arc of the radius R and the angle 2a (Fig. 3(a)), then Lq — 2RsiTia, L = 2Ra, 
and S ~^R^{2a — sin 2a), where the contact angle a is given by Eq. (|2^). Then, 

E = 2R{ea - Eq sin a) - ^jyR^ {2a - sin 2a) = (2a - sin 2a) (^eR - y^j^ i?^^ . (30) 

The height of the barrier is determined by the maximum energy at Rc — ec/^ajy: 

2 

Eb = {2a-sm2a)^-^. (31) 
2^ojy 

As one might expect, this barrier and consequently the relaxation rate do not depend on the film thickness. We stress 
that this estimation is valid only tor d > Rc. If d < Rc the energy barrier is obtained from Eq. (^^ by substituting 
R = d. This case of uniform current, however, leads to a thickness independent current density, and therefore cannot 
describe the experimental data. 



A. Surface pinning 

In this case, the whole current is confined to the surface layer of width A. It is apparent from Eq. (^ that for 
typical experimental fields (~ IT) A is smaller than the film thickness. This means that current flows mostly in a thin 
surface layer. Thus, all creep parameters, including the creep barrier, are governed by the total current ly, and not 
by the average current density Iy/2d. Then, apparently, the critical current density and the creep barrier are larger 
for thinner films, similar to the case of the collective- pinning effect mentioned above. Thus, also this scenario cannot 
explain the observed accelerated relaxations in the thinner films. 



B. Short-range bulk pinning 

Le t us c onsider the relaxation process for a critical state supported by the short-range pinning force discussed in 



Sec. [ID 2. The energy E of the vortex line is given by Eq. (ETp. The average critical current density corresponds to 



E = and is inversely proportional to the film thickness [see also Eq. (|2 

Ic ce . , 
^^^2d^2d$^^^""- (^'^ 

The energy barrier is given by the maximum energy atd — L + XftiL when the whole vortex line has left the potential 
wefi (Fig.4(a)): 

Eb = tana ^2c?esina — M^—ja^ , (33) 

where ja = Iy/2d is the average current density. If jc > ja > jc/2, then dE\j/dd < 0, i.e., the barrier is larger for 
thinner films. But, for ja < jc/2 the derivative dEt/dd > , and the barrier increases with the increase of the film 
thickness. Thus, under this condition {ja < jc/2) the magnetic relaxation rate is larger in the thinner samples. 



7 



The above analysis did not take into account the possibihty for dense defects. By "dense" we mean that the distance 
Ti from the neighbor potential well is less than dtana (see Fig. 4(b)). In this case the maximal energy (the barrier 
peak) is smaller than the barrier calculated in Eq. (|33|). Then the barrier energy is given by 

Eb = n l^e sin a - ^d^J^ (34) 

In this case dEh/dd < and the energy barrier for thinner films is always larger. Therefore one can see faster 
relaxation in thinner films only if the films are so thin that d < r^/ tana and the energy barrier is given by Eq. (|3^). 
From the experimental results shown below we infer that the average distance between effective defects > 1000 A 
in agreement with direct measurements using atomic force microscopy. 

To conclude, if the average current density in thin films becomes small enough compared to the original critical 
current density and if the films are thin enough, the relaxation at the same average persistent current is predicted to 
be faster for the thinner films. 



C. General case 



In the simplified picture of the critical-state relaxation outlined in the previous subsection, the total current was 
concentrated within a very thin layer of the width A. It was based on the assumption that the pinning force disappears 
when the vortex line leaves the small-size potential well, whereas inside the potential well the pinning force is very 
strong. As a result, outside the thin surface layers of the width A the vortex line consists of two straight segments 
(Figs. 3(b) and 4). In the general case, the distribution of the pinning force may be smoother and the shape of a 
vortex line is more complicated. In addition, interactions between the vortices may modify the barrier for fiux creep 
as well. However, the tendency must be the same: the current confined in a narrow surface layer drives the end of a 
vortex line away from the potential well to the regions where the pinning force is weaker and the vortex line is quite 
straight with the length proportional to thickness of the film if the latter is thin enough. Therefore, the barrier height 
for the vortex jump is smaller for smaller d. 

We also note that we do not consider an anisotropic case and limit our discussion to isotropic samples. The effect 
of anisotropy on the barrier height was considered in details in Ref. |l9|| In the presence of anisotropy the circular 
loop becomes elliptic and the vortex-line tension e must be replaced by some combination of vortex-line tensions for 
different crystal directions. These quantitative modifications are not essential for our qualitative analysis. 

Our scenario assumes that the current is concentrated near the film surfaces. In general, width of the current layer 
may vary from A to effective Campbell length Ac- One may then expect a non-monotonous thickness dependence 
when Ac is comparable with d. As we see, the Campbell length is an important quantity in determining whether 
current density inhomogeneity must be taken into account or not (in the absence of the surface pinning). The length 
Ac can be estimated from the micro-wave experiments: according to Golosovskii et al. |2l| Ac — 1000-\/i? A, where 
the field H is measured in Tesla. For H ~ 0.2 T this results in Ac ~ 450 A or 2Ac ~ 900 A, which has to be 
compared with the film thickness. 



IV. COMPARISON WITH THE EXPERIMENT 



A decrease of the measured current density with an increase of the film thickness is reported in numerous experimen- 
tal works P,p^ . This is consistent with the predictions given above for either surface or/and bulk pinning. Both 
pinning mechanisms predict similar 1/d dependence of j and it is, therefore, impossible to distinguish between surface 
and bulk pinning in this type of measurements. Only the additional information from the thickness dependence of 
the relaxation rate allows the drawing of some conclusions about the pinning mechanisms. 

Magnetic relaxation measurements in films of different thickness are discussed in detail in ^1 
from the data reported there we demonstrate an agreement of these data with our theory. 

Measurements were conducted on four 5x5 mm^ YBa2Cu307-s films of thickness 2d = 800, 1000, 2000 and 
3000 A, prepared by the laser ablation technique on SrTiOz substrates jl^. All samples had Tc « 89 K. The 
morphology of the samples was examined by atomic-force microscopy (AFM) technique and was found to be similar: 
the average grain size (1 — 50) x 10^ A and intergrain distance 50 A (For typical AFM picture of our samples, see Fig. 
1(c) in |12(| ). The magnetic moment was measured as a function of field, temperature and time, using a Quantum 
Design SQUID magnetometer. 



12 1 . Using excerpts 
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The average persistent current density was extracted from the magnetic hysteresis loops using the Bean model 
adapted for our case: ja [A/cm'^'^ = 30M/da^, where M [emu] is the irreversible magnetic moment, d [cm] is a half of 
the film thickness and a = 0.5 cm is the lateral dimension. Fig. 5 shows the persistent current density j a± T — b K 
as a function of the applied magnetic field H . Apparently, j is larger in thinner films. The same trend is found 
at all temperatures. These observations are in good agreement with Eqs. ( p^ and (p6|). We note, however, that 
since the value of js is not known, we cannot point out the dominance of pure surface, pure bulk or a mixed type of 
pinning. On the other hand it is unlikely that the observed thickness dependence is due to changes in the density of 
pinning centers with thickness, since the films' morphology is similar for all of our samples. This is further indirectly 
confirmed by the relaxation measurements. The decrease of current density due to increase of a mean grain size in 
thicker films would simultaneously result in faster relaxation, contrary to our observations. 

Fig. 6 shows typical relaxation curves at H = 0.2 T (ramped down from 1 Tesla) measured in films of different 
thickness. The interesting and unexpected feature is that curves cross, i. e., the relaxation is faster in thinner films. 
This is further illustrated in Fig. 7 where j vs. d is plotted at different times. At the beginning of the relaxation process, 
the average current density in the thinner films is larger. However, in the thinner films, the current density decreases 
much faster than in the thicker ones; as a result ja exhibits a non-monotonous dependence on thickness at later times, 
as s how n in Fig. 7. The faster relaxation in t hinne r films is in qualitative agreement with our results, discussed in 
Sec. in, in particular in subsections [II B and QIC. There, we find that such acceleration of the relaxation in thinner 
films may be understood only if we consider inhomogcncous bulk current density. In reality, it is very probable that 
both surface and bulk pinning mechanisms lead to inhomogcncous current density with a characteristic length scale 
in between the short (surface pinning) length A and the larger Campbell length. 



V. SUMMARY AND CONCLUSIONS 



Based on the two mode electrostatics approach we built a consistent theory of the critical state in thin typc-II 
superconducting films throughout the film thickness. We show that, irrespective of the pinning mechanism, current 
density is always larger near the surface, and decays over a characteristic length scale, which is in between A (of order 
of the intervortex distance) and the Campbell length A^. The length scale A is determined by the (finite) vortex 
tension and by the boundary conditions which force vortices to be perpendicular to the surface of superconductor, 
whereas the Campbell length Ac is determined by bulk pinning potential. 

Following this novel physical picture we conclude that: 

• Current density and magnetic induction in thin films in perpendicular field are highly inhomogcncous throughout 
the film thickness. Surface pinning significantly enhances these inhomogeneities. 

• Average current density decreases with the increase of film thickness approximately as 1/d. 

• Magnetic relaxation is slower in thinner films in the following cases: (1) In the absence of bulk pinning, i.e., 
only surface pinning is effective. (2) In the presence of bulk pinning, provided that the ratio between thickness 
and distance between neighboring defects is above a certain threshold {d/a ~ 1). 

• Magnetic relaxation is faster in thinner films only if bulk pinning is effective and the ratio d/a is below this 
threshold. 

In the experimental data presented here the measured average current ja decreases with the increase of film thickness 
as predicted, and the relaxation rate is larger for the thinner films, suggesting that d/a ^ 1, and the effective distance 
between defects > 1000 A. 
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Figure captions 

Fig.l Current density distribution vs. normalized depth z/d for the indicated d/\c ratios. The current distribution 
becomes more inhomogeneous as the ratio d/Xc increases. 

Fig. 2 Vortex energy (per unit length) in the vicinity of the pinning center of radius r^. 

Fig. 3 (a)Vortex depinning by a uniform current, (b) Simple scenario of vortex depinning by nonhomogeneous 
current flowing in a layer A. 

Fig. 4 (a) Barrier maximum configuration in the case of a dilute pinning centers, (b) Barrier maximum configuration 
in the case of the dense deffects. 

Fig. 5 Average persistent current density ja as a function of magnetic field at T = 5 ii" for films of difi'erent thickness. 
Fig. 6 Time evolution of the average persistent current density ja ai T = Tb K for films of different thickness. 
Fig. 7 Thickness dependence of the average persistent current density ja at T = 75 K taken at difi'erent times. 
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Figure 1 Prozorov et al. 




Figure 2 Prozorov et al. 




Figure 3(a) Prozorov et al. 




Figure 3(b) Prozorov et al 




Figure 4(a) Prozorov et al. 
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Figure 5 Prozorov et al. 
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Figure 6 Prozorov et al. 
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Figure 7 Prozorov et al. 



